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Abstract 



Recently, a novel class of models for inflation has been found in which the 
inflationary dynamics is driven solely by (non-canonical) kinetic terms rather 
than by potential terms. As an obvious extension, we show that a scalar field 
with non-canonical kinetic terms alone behaves like an energy component 
which is time- varying and has negative pressure presently, i.e. quintessence. 
We present a model which has a constant equation of state, that is, a "kinetic" 
counterpart of the Ratra- Peebles model of a quintessence field with a potential 
term. We make clear the structure of the phase plane and show that the 
quintessential solution is a late-time attractor. We also give a model for the 
"phantom" component which has an equation of state with w = p/p < —1. 
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I. INTRODUCTION 



Recent indirect or direct observations suggest that the Universe is currently dominated by 
an energy component with negative pressure |l]-f|]. One possibility for such a component is 
the cosmo logical constant. Another possibility is dynamical vacuum energy or quintessence, 
a temporary decreasing and spatially inhomogeneous component with negative pressure 
|5Hl8|. Only recently, a more radical candidate (called "phantom" component) has been 
proposed which is "growing" in time |TP| . 



Two problems arise from such a vacuum energy. The first is the fine-tuning problem: 
the vacuum energy density of order ~ 10~ 47 GeV 4 requires the introduction of a new mass 
scale about 14 orders of magnitude smaller than the electroweak scale. The second is the 
coincidence problem: the conditions in the early universe have to be set very carefully in 
order for the energy density of the vacuum and that of the matter to be comparable today. 
These problems are degenerate for the cosmo logical constant, however, they are separated 
in quintessence. A class of quintessence can avoid the coincidence problem by means of 
the attractor solution ||7|JT3|1. It is shown that the quintessence field approaches a common 
evolutionary track for a very wide range of initial conditions, so that the cosmology is 
extremely insensitive to the initial conditions. 

Usually the quintessence field is modeled by a scalar field with a canonical kinetic term 
and a potential term. However, we show that a scalar field with solely kinetic terms can 
(even without potential terms), albeit they are non-canonical, mimic such a (canonical) 
quintessence field. Our model is a natural extension of the kinetically driven inflation model 
proposed recently f20fl . A mechanism is proposed by which a dilaton remains massless |2l 



Assuming universality of the dilaton coupling functions, it has been shown that the dilaton 
evolves cosmologically towards values where it decouples from matter (so called "Least 
Coupling Principle" ) |2T] . 

Unlike the usual potentially driven quintessence model which automatically satisfies the 
weak energy condition, allowing non-canonical kinetic terms enables us to model the missing 
energy component which violates even the weak energy condition (so called "phantom field" 
19| ) . Only recently, Caldwell draw our attention to consider more general equation of state 



with w = p/ p < — 1 As he noted, since such a "phantom" equation-of-state cannot be 
achieved with a canonical Lagrangian and Einstein gravity, considering the phantom field 
requires some extension: either (i) non-canonical Lagrangian or (ii) non-Einstein gravity (or 
both). The latter possibility seems unlikely since the deviation from general relativity at 
the present time is strongly constrained by the solar system experiments [J22|. Our attempt 



is minimal one: non-canonical kinetic terms without a potential term. We intend to develop 
a more general study by including potential terms as well in the near future. 

The organization of the paper is as follows. In Sec.|I|, we present our model. In Sec. pi 



we start to show the existence of a scaling solution with a constant equation of state, 
— 1 < w < 0. Then we show that the scaling solution is a late-time attractor by means of 
linear and numerical analyses. In Sec.[IV|, we give a model which has a scaling solution with 
a constant equation of state of w < —1, and show that the scaling solution is a late-time 
attractor. In Sec.|V|, we make a comment on the possibility of reconstructing Lagrangian 
through observational data. Sec.|Vl| is devoted to summary. 
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II. BASICS 



We consider the following action of a single scalar field minimally coupled with gravity: 

S = Jd 4 x^g[-^R + P(^^) + S B , (2.1) 

where k 2 = 8ttG and Sb denotes the action of the background matter and/or radiation. 
Following |2(| , for simplicity, we only consider Lagrangians which depend only on the scalar 
field and its derivative squared^ 

x = -lyvv^. (2.2) 

For the general action of the scalar field ( |2.1| ), the field equations are given by 

R,u - \g,uR = k 2 ( ^^ X) V^V^ + p(0, X)g, u + T*\ (2.3) 

where T^ u denotes the energy-momentum tensor of the background. Eq.( [2.3|) shows that 
p((f),X) in the action ( |2.1| ) actually corresponds to the "pressure", p^, of the scalar field 
|20| , while the energy density, p^, is given by p^ = IXdpjdX — p. Hence the extrema of 
p(4>, X) with respect to X correspond to the same equation of state as that of a cosmological 
constant: p<p+p<p = IXdpjdX = 0. 

We assume that the universe is described by a flat homogeneous and isotropic universe 
model with the scale factor a. The time coordinate is so normalized that a = 1 at present. 
The field equations are then 

- = — -(pB + 3p B + p + 3p ), (2.5) 
a 6 

p B = -3H(p B + p B ) =: -3H(1 + w B )p B , (2.6) 

where ps and ps are the energy density and the pressure of the background matter and/or 
radiation, respectively. 

Since we only consider kinetic terms, we must impose that the function p((j), X) vanishes 
when X — > 0.0 Near X = 0, a generic Lagrangian may be expanded as 

p(<f>, X) = K{<P)X + L{<P)X 2 + ■■-. (2.8) 



x We use the metric signature ( — h ++). 

2 This amounts to assuming some resolution of the cosmological constant problem. The situation 
is the same as assuming the minimum of the potential energy is zero in the (canonical) quintessence 
field with a potential term. 
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III. POWER-LAW KINETIC QUINTESSENCE 



To see the effect of non-canonical kinetic terms in a concrete matter, in this section, we 
shall concentrate on the simplest Lagrangian containing only <p 2 and 4 terms, namely 

p(0,X)=K(0)X + L(0)X 2 . (3.1) 

In order to realize a model with negative pressure, K and/or L should be negative (note 
that X > 0). However, for the positivity of p^ for large X, we assume that L is always 
positive. We thus consider the case of K < 0. By redefining the scalar field and working 
with new field variable such that 

we rewrite Eq.( |3.1|) as 

p(0,X) = /(0)(-X + X 2 ), (3.3) 

where = new , X = X new = (L/\K\)X i d and f(<f>) = K 2 ((j) \ d ) / L(<p i d ). We may regard 
Eq.( |3.3|) as the basic Lagrangian. Then the pressure, pq, and energy density, pg, of the 
quintessence is given by 

p Q = f(<P)(-X + X 2 ), (3.4) 
p Q = 1X% -p = f(cf>)(-X + 3X 2 ). (3.5) 



A. scaling solution 

We look for scaling solutions which keep wq = Pq/pq constant. Then from Eqs.( [3.4|) 
and ( ]3.5| ), X is also found to be constant: 

X = l=f2.. (3.6) 
During the matter or radiation dominated epoch (p# ^ Pq), Eq. (|2.7|) becomes 

^ = -t7T^3( 1 +^W (3-7) 



Substituting Eqs. (|3.5| ) and ( |3.6| ) into the above equation, we thus obtain 

f{<f>) oc (0 - ^)-2(i+-q)/(i+» s ) ; (3. 8 ) 

where 0* is a constant. For simplicity, we henceforth choose <fi* = 0. 

To summarize, for the scalar field model with the constant equation of state, wq, during 
the matter or radiation dominated epoch, the function f((j)) should take the form of Eq. (|3.8| ). 
Conversely, if the function /(</>) is given by 



4 



f(4>) « 0" 



(3.9) 



then there exists a scaling solution such that the equation of state is characterized by 

1 + w B )a 



1. 



(3.10) 



Hence if we require that wq < during the matter dominated epoch, then the exponent a 
should satisfy 



a < 2. 



(3.11) 



Note that the weak energy condition (wq > — 1) can be violated if a < 0. The stability 
against perturbations is signified by the "speed of sound" defined by f20|| 



Pq,x = P,x 
Pq,x P,x + 2Xp >xx 



(3.12) 



For the quartic model Eq.( [3.4|) , one finds c 2 s = (1+Wq)/(5 — 3wq), and the model is unstable 
for perturbations on all length scales if the weak energy condition is violated: wq < — 1. 
One may wonder this is always the case. However, this is not so. In fact, we can consider a 
more general Lagrangian of the form |2^ 



p(<f>,X) = f(<f>)g(X), 



(3.13) 



where g(X) is an arbitrary function of X. One can show the equation of motion Eq.( p.7| ) 
has a solution X = const, with the function f(<ft) being the same form as Eq^^) if X 
satisfies 



ding 



dlnX 



= 1 +wq 

X=X Wq 



(3.14) 



Note that Eq. (|3.12|) involves the second derivative g t xx(X ). Therefore, one can always 
arrange it so that c 2 s > for stability. This may be a concrete realization of the "phantom" 



field (or "growing lambda") using the non-canonical Lagrangian ||19|| . An example of such 
a phantom field will be considered in Sec.|V[ In this section, we shall limit ourselves to the 
case of < a < 2 so that — 1 < wq < 0. 



B. attractor structure 

In the matter or radiation dominated universe, there exists the scaling solution for the 
kinetic Lagrangian. Here, we show that this solution is an attractor of the equation of 
motion for the scalar field. When ps ^> Pq, Eq.( |2.7| ) becomes 



) + TTrj , (1 - 2 )0 + f f (2 - 30 2 ) 2 = 0. (3.15) 
t(l + w B ) 4/ ^ / 



Since this equation has reflection symmetry <fi <-» —<p, we mainly consider the case of <fi > 0. 
Following |7|j25l, we make the change of variables as 



5 



t = \nt, 
u = 



(3.16) 
(3.17) 



where the scaling solution, S , is given by 



2(1 -w Q ) 
\j 1 - 3w Q 

With these changes, Eq. ( p.l5| ) becomes 



t =: £ 8 t. 



u = v , 
i/ = — f + 

X 



1 



i - ?>e s (v + uf 

9 n 

- {-( V + U )+ e s (v + u) 3 } - ^{-2(v + uf + 3g(v + uf} 



W B 



Au 



(3.18) 



(3.19) 



(3.20) 



where ' denotes the derivative with respect to r. Then one can find three critical points; 
(w, v) = (0,0), (1,0) and (—1,0). The (—1,0) critical point corresponds to the scaling 
solution with negative amplitude, while the (0, 0) critical point is a trivial one in which 
X = 0. 

In order to study the stability near the critical points, we use linear analysis |7|,|25 



Perturbing about the scaling solution (u, v ) = (1 + 5u,0 + 5v) and keeping only the terms 
linear in 5u and 5v, Eqs. (|3.19| ) and ( |3.20|) become 



5u' = 5v, 

5v' = —Sv + - — Q (Su + 5v) 



l + w B "~~ ~' 2(1 + w B ) I -Si 
Then the eigenvalues of small perturbations are given by 



(3.21) 
(3.22) 



A: 



2w Q - 1 - w B ± ^(2w Q - 1 - w B ) 2 + 8(1 + w B ){w Q - (1 + w Q )/(5 - 3w Q )} 

2(1 + w B ) ' 



(3.23) 



The necessary and sufficient condition for stability is that the real part of the eigenvalues 
be negative. Note that because — 1 < wq < 0, the second term under the square root is 
always negative. Then the condition for the stability is just 2wq — 1 — w B < 0. Hence 
the (1, 0) critical point corresponding to the scaling solution is stable. In a similar way, the 
eigenvalues of small perturbations near the trivial solution is given by 



A: 



2w Q - 1 - w B ± J(2w Q - 1 - w B ) 2 + 4(1 + w B )(w Q - 1) 



2(1 +w B ) 

Then the (0, 0) critical point is also stable. 



(3.24) 
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C. numerical analysis 



We have shown that the scaling solution is stable for small perturbations. In order to 
analyze the phase plane, we solve Eqs. (|3.19|) and (|3.20|) numerically. The phase plane is 



shown in FIG. |I]for the case of a = 1 and wb = 0. This figure shows the attractor structure 
of the scaling solutions and the trivial solution: there are trajectories which converge on 
these solutions asymptotically. The boundaries between phase flow, lines (1) and (2) in 
FIG. [I], correspond to the lines where Eq. ( |3.20| ) is singular: 



- U ±J±^-. (3.25) 
V6(1-^q) 1 ' 



Eq.( |3.20|) is also singular on u = except for 




The trajectories in the region v < —u + y(l — 3wq)/6(1 — wq) and u > approach and 
pass through the point (u, v) = (0,t>_), then converge on the scaling solution with nega- 
tive amplitude. From Eqs.( |3.12"|) , the requirement for stability against perturbations is not 
satisfied in the region 

1 — 3wn , N 9 1 — 3u>o ,„ 

Q <(u + vf<— (3.27) 



6(1 -w Q ) 2(1 -w q ) 

These regions correspond to the shaded ones in FIG. [[]. 

The cosmological evolution of the scalar field is obtained by solving Eqs. Q2.5p and ( |2.7| ) 
numerically. The initial conditions for the scale factor are so chosen that the Friedmann 



equation (|2.4j) is satisfied. Those for the scalar field are chosen in the region u > and v > 



—u + J(l — 3wq)/2(1 — wq) so that c 2 s > initially. In FIG. ||, we show the time evolution 
of energy densities of radiation, matter and quintessence field for various initial conditions. 
The present density parameter of the i-th component fli is defined by i\o = ^ViY^-^o)- 
We choose = 0.25 and set a = 1. The figure shows that for a very wide range of initial 
conditions, the energy density of the quintessence field converges on a common evolutionary 
track. 



D. mass scale 

In the Ratra-Peebles model of the quintessence field driven by a potential term, one 
can choose a parameter with mass dimension in the potential term to be a typical particle 
physics scale [0 . We introduce a parameter with mass dimension in the kinetic Lagrangian 



as 



f(<P) = -r^, (3-28) 
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then we fix this parameter by requiring that the scalar field is beginning to dominate the 
energy density of the universe today (p° rit is the present critical density and O is the present 
value of the scalar field) 

P°Q = /(0o)(-X + 3X 2 )^p° rit , (3.29) 
and that the scalar field has already reached the attractor solution 



L 
f 



-j- ~ H . (3.30) 



The last condition is obtained from Eq.( |3.15| ). These conditions fix the mass parameter and 
the present value of the scalar field as 

M ~ io( 43 «- 48 )/( 4 -°) [GeV], (3.31) 
0o ~ 10 43 [GeV" 1 ], (3.32) 

respectively. M could be larger than TeV scale when a ^ 1.3. 



IV. PHANTOM FIELD 



In this section, we construct a model of a power-law phantom field which is recently 
proposed as an alternative to the dark energy component |19[. In fact, in Ref. P9| , a 
scalar field with a kinetic term of inverted sign is introduced as a toy model of such a 
component: the Lagrangian density is given by Cp = +|V^0V M — V(<j>). The concern is 



the tachyonic instability for 0. However, as shown in [p!9| , as long as is negative, such 



an instability is not developed. For example, for a constant equation of state, one can show 



V,4>(t> — |(1 — w p)[H — t;H 2 (1 + wp)}, thus it is negative as long as —2 < wp < - 
model will not suffer from such a restriction even for a constant equation of state. 



-1. Our 



A. scaling solution 

In the previous section, we have mentioned that a general Lagrangian of the form 
Eq.( |3.13|) , p(4>,X) = f(<p)g(X), with the function f((f>) of the form Eq. fl3.9|) has scaling 
solutions with X =constant. From Eq. (|3.10|) , the equation of state is characterized by 

W P EE g£ = (1 + ™* )a - 1, (4.1) 

Pp 2 

then the phantom field which violates the weak energy condition, wp < —1, corresponds to 
a < 0. 

We determine the function g(X) which has a phantom solution. We impose the following 
conditions on g(X) other than g(X) — > for X — * 0(vacuum triviality) and g(X) > for 
X — > oo(positivity of the energy density for large X): (i) the positivity of the energy density 
Pp = f(2Xg : x — g) > because we are considering the missing energy component; (ii) the 



S 



violation of the weak energy condition, pp + pp = 2fXg t x < 0; (iii) the stability against 
perturbations, c 2 = g,x/{g,x + 2Xg, X x) > 0. 

In terms of (so that X = 2 /2) these conditions are rewritten as: (i') g — g > 0, 

which is geometrically interpreted as the intersection of the tangent to the curve g = g(<p) 
at the point with the vertical axis is negative; (ii') <fig ^ < 0; (iii') gx/{4>gxx) > 0. Thus, 
for > 0, g{4>) should satisfy g ^ < and g ^ < 0. Therefore, g(<fi) must have at least three 
extrema for > and hence be an eighth-order function of 0. We parameterize g{<p) as[]: 



n 



2 -6 2 )(0 2 -c 2 )(0 2 -d 2 )d0 



=: l0 8 - l -A^ + - 1C0 2 , (4.2) 

where b, c and d are constants corresponding to the extrema of g(<p) as shown in FIG. |3|. A 
sketch of g(<p) which satisfies the above conditions (i'), (ii") and (iii') is shown there. In the 
shaded regions, c 2 < so that the solution is unstable for perturbations on all length scale. 
From the above equation, the scaling solutions satisfy 

C0 2 /2-E0 4 /4 + A0 6 /6-0 8 /8 

Wp = : : : ; . 4.0 I 

C(p 2 /2 - 350 4 / 4 + 5A0 6 /6 - 70 s /8 
Then there are six scaling solutions at most. 



B. attractor structure 

We show that the phantom solution is an attractor of the equation of motion for the 
scalar field by means of the linear analysis employed in the previous section. From Eq.( |4.2j ), 
the equation of motion (|2.7D becomes 



0(C - 350 2 + 5A0 4 - 70 6 ) + — — -(C - 50 2 + A0 4 - 6 

t(l + w B ) 

f fl/~<i2 ^ ; 4 5 - 6 7- g 



+j [ifr - ^ + -a^ - -0« ) = o, (i.i ) 



when p B ^> p P . We make the change of variables as (|3.16| ) and 



u = — =:— (4.5) 
(p P £ P t 



where £p is determined by one of the solutions of Eq. (|4.3|) corresponding to the phantom 
solution. With these changes, Eq.(|4.4|) becomes 



3 We take negative sign for the coefficient of 2 as a minimal extension of the quartic model 
Eq. (|3.3|) . If we take positive sign, then the polynomial should be at least of tenth-order. 
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u 



v, 

-v + 



(4.6) 



1 



x 



C - 3B&(v + uf + 5A£j,(v + uf - 7$>{v + uf 
1 {-C{v + u) + B&(v + uf - A&(v + uf + Z%(v + uf} 



1 + wb 



-{- l -C{v + uf + \b&{v + uf - lA&(v 
u 2 4 o 



u 



huv + uf} 



(4.7) 



where ' denotes the derivative with respect to r. Then one can find seven critical points at 
most; (u, v) = (0, 0), (±1, 0) and others. 

Perturbing about the phantom solution (u,v) = (1 + 5u, + Sv) and keeping only the 
terms linear in 5u and 5v, Eqs.(O) and ( |4.7| ) become 



5u' = 6v, 

Sv' = —6v H Wp (Su + 8v) 

1 + wb 



2(1 + w P ) 
1 + w B 



S5u, 



where 



S 



C/2 - 35^/4 + 5Ae P /Q - 7£, 



C-3S^ + 5^-7^ 



& 9„ 



Then the eigenvalues of small perturbations are given by 



(4.8) 
(4.9) 



(4.10) 
(4.11) 



A 



2w Q - 1 - w B ± ^(2w P - 1 - w B f + 8(1 + w B ){w P - (1 + w P )S} 
" 2(1 + w B ) 



(4.12) 



Since we impose pp > and g < on the phantom solution, S 1 < 0. Then the second 
term under the square root is negative. Therefore the (1, 0) critical point corresponding to 
the phantom solution is stable. In a similar way, the eigenvalues of small perturbations near 
the trivial solution (0,0) is given by Eq.( |3.24| ) but wq is replaced with wp since the phase 
structure around the trivial solution is determined by the lowest kinetic term. The trivial 
solution is also stable. 



C. numerical analysis 

In FIG. U the phase plane for the case of a = — 1 and wp = is shown by solving 
Eqs.fl4~6l) and fl47| ) numerically. We choose b 2 = 1/2, c 2 = 1 and d 2 = 2(or A = 7/2, B = 
7/2, C = 1). There are one trivial solution and six scaling solutions. The shaded regions 
correspond to those in FIG. |3|, c 2 < 0. Eq. ( f4.7|) is singular on the dashed lines except 
for the points represented by cross. The boundaries between phase flow correspond to 
these singular lines. The trajectories in the region above line (1) and u > converge on 
the scaling solutions. In the region between line (1) and (2), they converge on the trivial 
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solution. In the region below line (2) and u > 0, trajectory approaches a point on the 
boundary asymptotically. 

The cosmological evolution of the scalar field is obtained by solving Eqs.Q2.5p and ( |2.7|) 
numerically. In FIG. [5], we show the time evolution of energy densities of radiation, matter 
and phantom field for various initial conditions. We choose Qm,o = 0.25 and set a = — 1. 
Since the scaling solution is growing in time and bounded by the instability regions (in 
which w > —1) as shown in FIG. the initial energy density should be chosen to be 
below the present energy scalefj Hence, there remains a severe fine-tuning problem of 
the initial conditions. We expect that including potential terms or introducing additional 
fields may allow initially decaying solutions. We note, however, that the condition (ii), 
p + p = 2fXg t x < 0, remains unchanged even if we include potential terms. 



V. RECONSTRUCTING P(<f>,X) 

It has been shown that the effective potential of a scalar field with a canonical kinetic term 
can be determined by using the magnitude-redshift relation of distant type la supernovae 
23| , |2"4|| . In a similar way, we consider the possibility of reconstructing p as a function of 



and X through observational data. The pressure and density of the quintessence or phantom 
component can be written in terms of the coordinate distance to redshift z, r(z) (which is 
related to the luminosity distance, (Il{z) = (1 + z)r(z)) as 

= X) = - f 2 - 2(1 + z) f] r,d f (5.1) 

(dr/dz) {dr/dz) 

* 2 P<f> = [2X^ - = 3 - 3H 2 n M , (l + z)\ (5.2) 

From Eq.( |5.1| ), we find that p(<fi,X) is written as a function of z. In order to reconstruct 
p(4>,X), we further need to rewrite (or X = 2 /2) as a function of z. One might expect 
that this may be done using Eq.( |5.2| ). Indeed, it is possible for the case of the quintessence 
field with a canonical kinetic term and a potential term, where p(<f),X) = X — V{(f) and 
p(4>, X) = X + V(<f>). By summing up Eqs.(|5.1|) and (|5.2|), X can be described as a function 
of z. Integrating X with respect to z enables us to rewrite as a function of z. Then, 
by subtracting Eq.( |5.1|) from Eq. (|5.2|) , V(4>) can be described as a function of z, that is, 
cf). However, in our case, p(4>, X) is an arbitrary function of both <p and X so that it is 
impossible to rewrite <fi or X as a function of z without identifing the combination of and 
Xinp(0,X).Q 

Of cource, once one specifies the functional form of p(<f), X), the reconstruction is possible 
like the case of the quintessence field with a canonical kinetic term and a potential term. For 
example, for the case of our model Eq. (|3.3|) , we can write f((p) and X in terms of observable 



*It is even possible to start from zero energy density with a fine-tuning of the initial conditions. 



5 We also note that the reconstruction method fails if the scalar field has a multi component or is 
non-minimally coupled 1 15. 16] to the curvature. 
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quantities in the following manner: Take the ratio of Eq.( |5.2|) to Eq.( |5.1|) . Then X is written 
as a function of z. Integrating X with respect to z, <fi can be written as a function of z. 
Once again using Eq. ([5TT|) together with the results obtained, we can rewrite / as a function 
of z. Hence we can reconstruct /, /(</>) = f((f>(z)). 

VI. SUMMARY 

We have shown that a scalar field with only non-canonical kinetic terms can, without 
the help of potential terms, behave like quintessence or phantom energy component. 

We have presented a kinetic counterpart of the Ratra-Peebles model and investigated 
the structure of the phase plane and shown that the quintessential solution is a late-time 
attractor. The phase area of the initial conditions for the scalar field which converge on 
the quintessential solution is smaller than that of the Ratra-Peebles model since there are 
a trivial solution as well as the scaling solutions in the kinetic Lagrangian. However, for 
very wide ranges of initial energy density, these energy components converge on common 
evolutionary tracks. 

We have also given a model of a power- law "phantom" field with arbitrary w of w < — 1 
which is stable against perturbation, and have shown that the phantom solution is a late- 
time attractor. In our model of the phantom field, the initial conditions have to be set 
carefully because there are a number of scaling solutions. Moreover, the initial energy 
density should be chosen to be below the present energy scale because the decaying lambda 
region corresponds to the unstable region. Therefore, there remains a fine-tuning problem 
of the initial conditions, which is no more severer than that of the cosmo logical constant. 
The situation may be improved by including potential terms or by introducing additional 
fields. 

The violation of the weak energy condition is required to construct wormholes |26) . Our 
model of a phantom field may be used to explore such a fascinating possibility. 
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FIGURES 




U = 0/0 s 

FIG. 1. Phase plane of the quintessence field in the matter dominated universe. Trajectories 
in the region between line (1) and (2) converge on the trivial solution (u, v) = (0, 0). In the region 
above line (1) and u > they converge on the scaling solution (1,0) directly, while in the region 
below line (2) and u > they approach and pass through the point (0, v-), then converge on 
another scaling solution (—1, 0). The speed of sound is imaginary in the shaded region. 
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6 4 2 

log 10 (l+z) 

FIG. 2. Energy densities of radiation, p/j, matter, pm, and quintessence field for various initial 
conditions, Pq , i = 1,2,3, against the redshift. The initial conditions for p^ correspond to 
that of the scaling solution in the radiation dominated universe. The initial energy density of 
the quintessence field which converges on the attractor solution spans more than 10 orders of 
magnitude. 
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FIG. 3. A sketch of the function g((j>). c 2 s < in the shaded regions. The phantom solution is 
one of the solutions of Eq.(fO|) at which <pg ^ < and g ^/<pg ^ > and the intersection of the 
tangent with the vertical axis is negative. 
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-1 1 (2) 

U = 0/0 p 

FIG. 4. Phase plane of the phantom field in the matter dominated universe. We set b 2 = 1/2, 
c 2 = 1 and d 2 = 2. The scaling solutions are represented by open circle. Eq.( [4.7D is singular on 
the dashed lines except for the points represented by cross. In the shaded regions corresponding to 
those in FIG|], c 2 < 0. Trajectories in the region between line (1) and (2) converge on the trivial 
solution. In the region above line (1) and u > 0, they converge on the scaling solutions. In the 
region below line (2) and u > 0, trajectory approaches a point on the boundary asymptotically. 
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6 4 2 

log 10 (l+z) 

FIG. 5. Energy densities of radiation, pR, matter, pm, and phantom field for various initial 
conditions, pp , i = 1,2,3, against the redshift. The initial conditions for p^ correspond to that 
of the phantom solution in the radiation dominated universe. The initial energy density of the 
phantom field which converges on the phantom solution spans more than 10 orders of magnitude. 
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